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C ! 1. INTRODUCTION 

^ , 

QQ ' Over the last twenty years, widely different conjectures have been proposed for the infrared behaviour 
. of the gluon propagator. Although it is a gauge dependent quantity, it can be discussed in a given 
gauge. Even within the same gauge, the proposals for the infrared dependence differ drastically [Q. We 
mainly summarize here the results that are given in the literature within the Landau gauge, since that 
' gauge is widely used in studies of Dyson Schwinger equations (DSE) as well as in lattice QCD. Early 
■ predictions were obtained by solving approximately the DSE. Mandelstam |^ obtained a solution of a 
, set of truncated DSE equations with an infrared behaviour of the form (g^)^^ for the gluon propagator. 

Such an infrared enhancement was shown, if obtained in any gauge, to lead to an area law for the 
. Wilson loop [H and thus to be sufficient for confinement. Infrared enhancement was assumed in various 
I phenomenological studies Q and corroborated by later studies of DSE with refined approximations |^ . 
A different perspective was taken by Gribov who showed that avoiding gauge copies one would 
obtain a gluon propagator which vanishes in the infrared in the Landau and Coulomb gauges, of the 
pi • form 
(D ' 
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D{q')-^^ ■ (1-1) 

, An infrared suppressed behaviour was advocated by Stingl 0, and recently by others as a possible 
^ ■ solution to DSE. Following a procedure similar to that by Gribov, Zwanziger |^ gave arguments to show 
that, on the lattice, for any finite spacing in the limit of infinite volume, D[q^ = 0) = 0. 

We will also consider in this work the parametrization deduced by Cornwall [ p^ using a resummation 
of Feynman graphs which leads to gauge-invariant amplitudes. The gluon propagator is obtained as a 
solution to this special set of DSE where the claim is that the only gauge dependence appears in the 
free part. Cornwall's solution, in addition to fulfilling the Ward identities, allows a dynamical mass 
generation. Thus this formulation has the additional attractive feature that the gluon mass vanishes in 
the ultraviolet as required perturbatively. Since the self-energy obtained by Cornwall is claimed to be 
gauge-independent, we will use his model to fit the propagator both in the Landau and in the Laplacian 
gauge. 

In contrast to all the approaches described above, lattice QCD provides a framework for the calculation 
of the gluon propagator starting directly from the QCD Lagrangian and can thus yield a conclusive 
result. Attempts to calculate the gluon propagator started more than ten years ago |ll|,|l^ on rather 
small lattices. These early results could be interpreted in terms of a massive scalar propagator, but 
confirmed the expectation that a Lehmann-Kallen representation is not applicable: positivity of the 
transfer matrix is lost after the non-local gauge fixing. Results on larger lattices were accounted for by 
assuming a positive anomalous dimension [|l3|. Recently, a detailed study of the gluon propagator on 
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very large lattices jlj] has been performed, which makes an impressive effort towards bringing under 
control errors due to the finite lattice spacing and to the finite lattice volume. However, up to now, all 
lattice studies have used a similar implementation of the Landau gauge on the lattice. Gauge-fixing is 
accomplished by using a local iterative procedure which identifies local stationarity, but in general fails 
to determine the global extremum. Which local extremum ("lattice Gribov copy") is selected depends 
on the starting condition. These lattice Gribov copies cannot be eliminated. In this situation, their 
effect has repeatedly been claimed to be small |T^]. As discussed in Section 3, we are not convinced by 
such claims. Results obtained so far will have a safer foundation if the effects of lattice Gribov copies 
are better understood. 

In this work we address the problem of Gribov copies. We use a different gauge condition, which 
produces a smooth gauge field like the Landau gauge, but which specifies the gauge uniquely: no 
ambiguity arises due the lattice gauge fixing procedure. This is accomplished by using the Laplacian 
gauge Ip^ . The motivation and implementation of this gauge are given in section 3. 

We calculate the gluon propagator in quenched QCD on lattices of sizes 8'*,I6'' and 16^ x 32 at 
(3 = 5.8 and 6.0, in an attempt to study its zero-temperature behaviour. Our procedure can be extended 
straightforwardly to finite temperature where the infrared behaviour of the propagator yields the chromo- 
electric and chromo-magnetic screening masses. The results that we obtain, within the Laplacian gauge, 
show the same ultraviolet behaviour as in Landau gauge. However, there are significant modifications 
in the infrared. In particular we find that the zero-momentum propagator is finite, obeys scaling, and 
becomes volume independent for large enough volumes. It should not however be used as a definition 
of the gluon mass, since the zero-momentum limit of the propagator is gauge dependent. It is simply 
a measure of the susceptibility of the gauge-fixed field in the Laplacian gauge. A quantity which 
instead can be shown to be gauge independent to all orders in perturbation theory is the pole mass of the 
transverse part D(q^) of the propagator To determine this pole if it exists at all, an extrapolation 
to negative is necessary. We compare the inverse propagator D~^{q^) in the Laplacian and the 
Landau gauges. Using a variety of extrapolation ansatze, in particular a fit to Cornwall's model 
which describes the momentum dependence of our results rather well, we find that data in the Laplacian 
gauge give support for the existence of a pole at a mass of ~ 640(140) MeV. Data at smaller momenta 
are needed to consolidate this result. 

Section II introduces our notation; Section HI motivates and describes our choice of the Laplacian 
gauge; Section IV presents our results. They are summarized in Section V. 



2. DEFINITION OF THE GLUON PROPAGATOR 



The gluon propagator in the continuum is given by 

K'li'l) = d^a:(0|T[A;i(x)At(0)]|0)e*«- (2.2) 

This tensor can be decomposed into a transverse and a longitudinal part: 

<(g) = (^Sf,. - ^-if-^ S^' D{q-) + ^-^6^"^ (2.3) 

For a covariant gauge F{q^) reduces to a constant and corresponds to the gauge fixing parameter ^ 
which in the Landau gauge is zero. Since we want to make a comparison with the recent results 
obtained in the Landau gauge, we study the transverse scalar function D{q^) which can be extracted 
from V-^M- 

^^'^')-\{y.\i:k»]-\^ • (2-4) 

F[q^) is determined by projecting the longitudinal part of ^"^(q) using the symmetric tensor q^q'^ . 
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On the lattice the dimensionless gluon field can be defined by 

A^{x + A/2) = ^{ [t^p(^) - Ki^)] - ^Tr [U^{x) - C/t(x)] | + 0{a^) (2.5) 

where a is the lattice spacing. One may consider different definitions for the gluon field A^, accurate to 
higher order in a. It has been found [l9| ] that these different definitions give rise to modifications that 
can be absorbed in the multiplicative field renormalization constant. 

The gluon propagator in momentum space is constructed by taking the discrete Fourier transform of 
for each colour component 

= E ^-'"■^''^^'^^K^^ + A/2) (2.6) 

X 

where the discrete momentum q — (g^,/i — 1, ..,4) takes values 

9m = J^'^M' = -(^L^- 1)>->(^L^) (2-7) 

and the momentum-space gluon propagator D'^^^^{q) is defined by 

Y 6{q-q')D';i{q)^mQ)Al{~q')) • (2.8) 

with V the lattice volume. In the ultraviolet the gluon propagator is expected to behave like 1/g^ . Since 
on the lattice the free massless propagator behaves as 

Diq) = 2, (2.9) 

to reduce errors due to the finite lattice spacing we take as our momentum variable the usual 

2 a,, a , , 

= -sin^ (2.10) 

To relate the bare lattice propagator to the renormalized continuum propagator -DR(g; fi) one needs 
the renormalization constant Z3(/i, a): 

a^Diqa) = Z3i^^,a)D^i{q■,|I) . (2.11) 

Imposing a renormalization condition such as 

DR{q)U2^^2 = ^ (2.12) 

at a renormalization scale fi allows a determination of Z^{fi,a). Connection to other continuum renor- 
malization schemes can then be made. 



3. GAUGE FIXING PROCEDURE 
3.1. Motivation 

The gluon propagator is normally considered in Landau gauge, d^A^{x) = Va;. On the lattice, this 
condition becomes: 

F{n) = ^Rc Tr {n{x)^U^,{x)n{x + /i)) maximum (3.13) 

The gauge-fixing functional F has many local maxima. To specify the gauge uniquely, the gauge condi- 
tion above refers to the global maximum. This defines the Fundamental Modular Region (FMR) Landau 
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gauge. In practice however, the gauge transformation is found by an iterative local maximization of 
F, which terminates when any local maximum has been reached. A different gauge condition is thus 
implemented, which one might call the random Landau gauge, and which depends on the details of the 
maximization procedure. 

It is commonly believed that the effect of choosing a local maximum of (I3I3I) rather than the global 
maximum is small, so that the "random" Landau gauge is a good approximation to the FMR Landau 
gauge. The following argument is often presented to support this view. A given gauge configuration 
is gauge-fixed n times, each time after performing a random gauge transformation; this procedure 
generates many gauge copies, each corresponding to the local maximum nearest to the random starting 
point along the gauge orbit. It is observed |0 that the difference between gluon propagators measured 
on copies corresponding to the largest and the smallest values of ( ^.13 ) is found to be statistically 
insignificant. A possible problem with this argument however is that the number n of gauge copies 
considered in such comparisons (typically 30 or less) is extremely small compared to the total number 
of local extrema of (3.13): for simple entropic reasons, all copies considered miss the global maximum 
by similar amounts, and no reliable information can be extracted about the gluon propagator in the 
global maximum configuration. It is therefore possible, and we believe quite likely, that the "random" 
Landau gauge and the FMR Landau gauge are significantly different. 

Further evidence for this situation has recently been provided in another gauge, the Direct Maximal 
Center (DMC) gauge Although the functional FjjMci^) to be maximized differs from (3.13), a 

similar approach of local iterative maximization is taken, leading to the "random" DMC gauge, with 
similar problems. In this case however, it is also possible to converge to a large value of Fdmc by 
starting from a Landau gauge copy ( "Landau" DMC gauge) . This value F]^ can then be compared with 
the values obtained from n random starting points. One may fit the maximum value among n copies, 
F{n), by a reasonable ansatz like a series in and extrapolate to n ^ 00. It turns out that the 
extrapolated value falls well below Fl, which is itself below the global maximum j2^. Furthermore, the 
properties of the gauge-fixed field are qualitatively different between the "random" and the "Landau" 
DMC gauges: the former confines after center projection, while the latter does not p6| . 

Since in Landau gauge as in DMC gauge, the number of local maxima is expected to grow exponentially 
with the lattice volume, we expect a similar situation in Landau gauge, leading to large differences 
between the "random" (local maximum) and the FMR (global maximum) gauges. One might argue 
that this is not a problem, and that the local maximization of (3.12) implements in the thermodynamic 
limit a well-defined, but stochastic gauge condition. The relationship between that gauge condition and 
its perturbative version 9^Ap(x) = is unclear however. Therefore, one should consider the possible 
effects of selecting a local rather than the global maximum of ( 3.13| ) with a great deal of caution. This 
is the motivation for our study of the gluon propagator in a wcU-dcfincd, unambiguous gauge. 



3.2. 5(7(3) Laplacian gauge fixing 

In ll^, Vink and Wiese proposed a simple method to fix the gauge unambiguously in SU{N). It 
uses N auxiliary Higgs fields, which are chosen as the N lowest-lying eigenvectors w'^*-' of the covariant 
Laplacian. Under a local gauge transformation fl{x), these eigenvectors transform covariantly: u^*) (x) — > 
fi(x)i;*^*'' (x). Therefore, the gauge can be fixed by requiring, at each space-time point x, {f2(a;)w(*^ (x), i — 
1,..,A^} to have some predefined orientation in color space. Specifically, each eigenvector v'^'^\x) has 
N complex color components, so that the N eigenvectors form a complex N hy N matrix M. Ref. 
projects this matrix onto SU{N) by polar decomposition: M = WP,W G U{N),P = (A/tM)!/^. 
The required gauge transformation is then ^{x) = e'^°'W\ where a = -^arg(det W^). ^{x) rotates M 
"parallel" to the identity l^r at each space-time point. The gauge is unambiguously defined, except for 
these gauge configurations where some of the N lowest eigenvalues are degenerate. Such configurations 
are genuine Gribov copies; they never occur in practice. This approach has been tested for SU{2) 
and f/(l) [ p3| and it was shown to reduce to the Landau gauge in the continuum limit aside from 
exceptional configurations (e.g. an instanton background). Here, we use a slightly modified procedure 
which requires only (N — 1) eigenvectors (2 for SU{3)), as follows p^ . 
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First, apply a gauge transformation ri^^^(x) which rotates v^^\x) to . Five real compo- 




nents of the rotated w*^^-' (x) must vanish, which specifies five constraints. Therefore fi*^^^ {x) is not fully 
specified, but has 8 — 5 = 3 degrees of freedom. Any satisfactory D,^^^ can be used. 



To completely fix the gauge, we use the second eigenvector v^^\ already rotated by fi^^^ to 

„,(2) 



/ (2) 



Three additional constraints are obtained by requiring to be rotated to | |^(2) p |^,(2) p | . This 



fixes the gauge completely and uniquely. 

Note that the second rotation is in an SU{2) subgroup, since it leaves v^^ untouched. This indicates 
how to generalize this construction to SU{N): the first rotation fixes (27V— 1) constraints, which leaves 
{N'^ — 1) — (27V— 1) — {{N — 1)^ — 1) degrees of freedom, forming a subgroup SU{N — 1). The next step 
reduces the gauge freedom to SU{N — 2), etc... down to SU{2). It is easily seen that, in this recursive 
procedure, the matrix M is reduced to upper triangular form (with real positive diagonal elements) by 
the rotation i^{x). This is why the TV*'* eigenvector needs not be computed: it is only transformed by 
a phase, which is separately determined by the requirement that Q{x) £ SU{N). Our procedure can 
thus be viewed as a QR decomposition of M. The gauge, which is globally well-defined (provided the N 
eigenvalues are distinct), may be ill-defined on a sub-manifold of points x where our recursive process 
breaks down. It can be seen that such local gauge defects occur at isolated points, where for SU{3), 

\J ^2^'' P + ba^' P = 0. The correlation of these points with instantons is studied in |24|. 

The Laplacian gauge so defined has the great virtue of being unambiguous. Hence it is the appropriate 
tool to address our concern about the effect of local extrema of the usual Landau gauge. It also has 
strong similarities with Landau gauge: it is smooth, Lorentz-symmetric, and gauge-fixes a pure gauge 
lattice configuration (gauge-transformed from U^{x) — 1 Va;,/i) back to [/ = 1. Nevertheless, it is a 
different gauge: its perturbative definition is under consideration [25|] ; it differs from Landau gauge most 
strongly where the magnitude of the eigenvectors |v'-^'^^(x)| becomes small. 



4. RESULTS 



Since most of the previous studies were performed in the Landau gauge, it is important to compare 
our Laplacian-gauge propagator with the Landau-gauge one. For this purpose, we have taken, for our 
analysis, lattice configurations available on the Gauge Connection database which had already 
been gauge fixed to Landau gauge with the usual local Over- Relaxation method ||2|]. These are 200 
configurations of a 16^ x 32 lattice, at /3 = 5.8 and 6.0 each. 

The transverse gluon propagator is shown in Fig. |l| for the two gauges at /3 = 6.0. As expected 
the ultraviolet behaviour is identical in the two gauges, whereas in the infrared, which is the region of 
interest, significant differences are visible. Since we use a different gauge, this should not come as a 
surprise. We show the usual quantity q^D{q^). The Laplacian propagator is clearly not as large as the 
Landau propagator at low momenta. 

The difference between Landau and Laplacian gauge can also be seen in the deviation of F{q^) from 
zero. Whereas in Landau gauge we find that 

rq'^'Dll « 1 (4.14) 

as expected, in the Laplacian gauge F{q^) is not small, and has a maximum at low momenta. The 
behaviour of F{q^) is shown in Fig. ^ for 8^ and 16^ x 32 lattices at (3 — 6.0. Since F{q^ = 0) can not 
be obtained by our projection, we only have one point, at the smallest momentum 27r/32 on the larger 
lattice, to ascertain that F{q^) really has a maxinmm and does not keep diverging as ^ 0. But since 
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the data are systematically higher for the smaller lattice than for the larger one, it seems unlikely that 
increasing the lattice size further would bring the infrared data up and remove the maximum. 

It is interesting to examine the volume dependence of the zero-momentum propagator ©(0) = 
2 X^a -^M/i ~ 0)- order to determine the transverse part of the propagator 

at zero-momentum, -D(O), one must subtract from V{q^) F{q^)/q^\q-2^Q, which we can only obtain as 
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FIG. 2. F{q^) as a function of |g| in lattice units after applying the cylindrical cut in momentum. 
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FIG. 3. The Renormalized zero-momentum propagator ©(0) = ■jX^/j^X/a ~ 0) versus volume in 

physical units. The dashed line is a fit to the form a exp ( — V/Vo) + c. 



/^From Fig. it can be seen that to extract the limit of F{q'^)/q^ as ^ reliably, one needs 
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more data in the infrared. Therefore, with the volumes at our disposal we can only examine the 
zero-momentum limit 2?(0) of the full propagator. In the Landau gauge, Zwanziger has argued that 
D{Q) — I?(0) should vanish in the infinite lattice volume limit A recent lattice study in SU{2) at 
finite temperature seems indeed to indicate such a behaviour. With the data on our present volumes 
the needed subtraction in the Laplacian gauge cannot be reliably performed, and thus we cannot extract 
£'(0). What we find is that I'(O), the zero-momentum propagator, is finite and volume independent 
for large enough volumes. The volume dependence and scaling of the rcnormalized zero-momentum 
propagator in physical units is displayed in Fig. ^ where we collected results from /3 = 5.8, 6.0, 6.2 and 
6.5. To obtain the rcnormalized propagator we impose the renormalization condition given in eq. ( ^.12 ) 



where we choose the renormalization point to be /i = a for /3 = 6.0 i.e. /i = 1.943 GeV. This 



determines Z3(/x, a^^g.o) ~ 2.312. We then use eq. (2.11) to find the ratio of the factors Z3 for different 



/3 values at the same physical momentum q, e.g. for /? — 5.8 a.t q = fi we have: 

Zsif^, a 13=5.8) _ °^^5.8^(g"/^=5.8) 
^3(^^,0^=6.0) a^=6.0^('?«/3=6.o) 

= 0.97(4) (4.15) 

In this way we obtain the renormalization factors at all (3 values. For f3 — 6.2 and P = 6.5 we find 
1.04(4) and 1.07(6) respectively as compared to the value at /? = 6.0. We also obtained consistent 
results by fitting our data in the ultraviolet regime using the asymptotic one-loop result for Dji ^ 
Z/g^(l/21n(g^/A^))~'^°, with djj = 13/22 as in the Landau gauge since in this regime the results in 
the Laplacian and Landau gauges are the same. As can be seen from Fig. ^ the rcnormalized I?(0) 
displays reasonable scaling, and appears quite volume- independent for volumes larger than 1/2 fni**. 
We find a value of 2?(0) = 5.020(16) GeV-^, or V{0)-^/^ = 445(3) MeV, corresponding to a length 
scale of ~ 0.5 fm. Since the zero-momentum propagator measures the susceptibility of the field, the 
length associated with it determines the domain over which the gluon field remains correlated in the 
Laplacian gauge. If the lattice dimensions become of the order of this characteristic length, then one 
expects finite size effects to become appreciable. This is indeed what is observed, as shown in Fig.||, 
with an approximate volume dependence of exp{—V/Vo) with V the lattice volume and Vq ~ I?(0)^. 

On the lattice, the Lorentz symmetry is only approximately restored. Lattice artifacts cause some 
dependence of D{q) on the orientation of the vector q rather than just on q^. To minimize these 
discretization effects, we filter our data by making a cylindrical cut in momentum along a reference 
direction ri = ^(1, 1,1,1), in the same manner as in Rcf. |l^ . Namely, we only consider momenta 
obeying the criterion \ A.q\ < 2t:/Ls, where Lg is the number of sites in the spatial direction, and Aq is the 
momentum transverse to h (Aq — q — q.n h). Using these filtered data which allow a direct comparison 
with [Q, we examine the various proposals discussed in the Introduction for the infrared behaviour of 
the propagator. We find that Gribov type parametrizations [^,0 as well as infrared enhancement of the 
type {q^y^ are excluded H]. The ansatz of Marenzoni et al. (l3), 

^('^'^ - WT^Tm^ ' (4.16) 

with a non-perturbative anomalous dimension a, gives a better description of the lattice data than the 
aforementioned parametrizations, but, as seen in Fig. ^, underestimates the peak of the propagator. 
On the other hand, Cornwall ||l0[| allows for a dynamically generated gluon mass which vanishes at 
large momentum in accord with perturbation theory. Using a special set of DSE referred to as a gauge 
invariant "pinch technique" , he obtains the following solution for the gluon propagator 



{q^ + M\q^))\n'i -j^-^ 



with 



,2, f In [(g2 + 4M2)/A21 1 ^ ^ 
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FIG. 4. The gluon propagator -D(q^) multiplied by on the 16^ x 32 lattice at /? : 
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i.O. The dashed-dotted 



line shows the fit to the model by Marenzoni et al. eq.(4.16), the solid line to Cornwall's model eq.(4.17) and 
the dashed line to model A of Ref. eg. (118). 



Cornwall's proposal provides a reasonable fit to the data over the whole momentum range (with 
X^/n.d.f — 2.5). The quality of this fit can be seen in Fig. |4[ For comparison we also fitted our data to 
the form suggested by Leinweber et al. [ [l^ where two terms were used, one to describe the ultraviolet 
behaviour of the form Z3tiv ~ j>. L{q'^,M), and one the infrared of the form Dis, ~ 1 / {q^ + Apy'^" . 



The exact form, referred to as model A, as taken from ref. [g^, is 



Die) 

L{q\M) 



AM 



,f2a 



-In [(q2+M2)(q- 



g2 + M2 



Hq^M) 

-13/22 



(4.18) 



This parametrization, which includes one more parameter than Cornwall's and is purely phenomenolog- 
ical, does fit the data best over the whole momentum range (with x^/^-d-i = 1-2). 
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FIG. 5. Scaling of the data at P — 5.8 and 6.0 on the 16'^ x 32 lattice. The solid curve is the best fit to both 
sets of data. 
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We address the question of scaling by comparing our results at /3 = 5.8 and (3 = 6.0 on the largest 
lattice. In the scaling regime the renormaUzed propagator I3r,((7; /i) is independent of the lattice spacing. 



Therefore, as in ref. we can use Eq. (2.11) to obtain the following expression for the ratio of 

unrenormalized lattice propagators at some physical momentum scale q 



-Pi(gQi) _ Z3{^J.,al)D^l{q■, ^)/al 
D2{qa2) Z3{fj,,a2)DTi{q;fi)/al 



(4.19) 



and where the labels 1, 2 refer to the data at /3 = 6.0 and (3 — 5.8 respectively. The scaling properties of 



the lattice gluon propagator can now be investigated using Eq. (4.19) by adjusting the ratios Z1/Z2 and 
ai/a2. In Fig. ^ we show the two sets of data lying on the best scaling curve. The shifts required along 
the horizontal and vertical axes determine the ratios of the wavefunction renormalization constants and 
of the lattice spacings. We find 



a/3=6.o/a/3=5.8 = 0.71 ± 0.02 



and 



^/3=6.o/^,3=5.8 - 1.07 ±0.075 



(4.20) 



with strongly correlated errors. The ratio of lattice spacings is in agreement with the value of 0.72(4) 
obtained from a detailed analysis of the static potential . The ratio of the Z-factors is within what 
is expected from perturbation theory, and in agreement with the value of 1.04(3) of Ref. ||lj]. In other 
words, scaling is very well satisfied for the Laplacian gauge, and performing the fits at P — 6.0 gives the 
behaviour of the gluon propagator in the physical regime. 

We focus now on the infrared behaviour of the transverse propagator. Figs. |^ and ^ show the inverse 
propagator as a function of q'^ in the two gauges. Two advantages of the Laplacian gauge become visible. 

First, the orientation of the momentum q has less effect than in Landau gauge: the data points at 
a given value of q^ show less scatter, and the cylindrical cut is not as essential as in Landau gauge in 
the infrared region. At a given lattice spacing, the Laplacian gauge approximates better the Lorentz 
symmetry of the continuum. This reduction of lattice artifacts is understandable since the gauge is 
fixed by considering the lowest-lying eigenvectors of the Laplacian, which are the least sensitive to UV- 
cutoff effects. In contrast. Landau gauge comes from the iteration of a completely local, UV-dominated 
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process. Better rotational symmetry allows for better accuracy, or for the same accuracy on coarser 
lattices. 

Second, the inverse propagator is closer to a linear function of in Laplacian gauge. If it were the 
propagator for a free boson, it would be described by a straight line since 1/D{q^) = Z~^{q^ + w?). 
Having curvature means that one has a momentum-dependent effective mass n(g^). In particular, the 
infrared mass 11(0) and the pole mass n(g^) such that {q^ + n(g^)) = are different. 
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FIG. 6. The inverse gluon propagator at low momentum in Landau gauge, at /3 = 6.0 on the 16^ x 32 lattice. 
The filled triangles and crosses show the data which arc kept and discarded by the cylindrical momentum 
cut respectively. Three extrapolations to negative are shown: quadratic and cubic polynomials in , and 
Cornwall's model. Note the instability of the pole D~^{q^) = with respect to the type of extrapolation chosen. 

The latter is of special interest, because of its gauge independence at least to all orders in perturbation 
theory. Finding a pole, i.e. a zero of the inverse propagator, requires the extrapolation of our data 
to negative ^. The less curvature in the inverse propagator in the infrared, the more reliable the 
extrapolation will be. 

Three types of extrapolation are displayed in the figures: quadratic and cubic polynomials in q^ , and 
our fit to Cornwall's model. The location of the pole, and even its existence, are affected by the choice of 
extrapolation in the Landau gauge. The coefficients ai, 02, 03 of the cubic polynomial extrapolation keep 
increasing, indicating poor stability. Essentially, no statement about a pole can be made in that gauge. 
Differentiating between a cubic fit (which gives a pole) and a Cornwall-type fit (which doesn't) will 
require extremely accurate data on large lattices. Extracting the pole is also difficult in the Laplacian 
gauge but at least one finds a pole with all the ansatze that we tried. Given the convexity of the data, a 
lower bound is provided by a linear fit near q^ = 0, which defines the (gauge-dependent) infrared mass. 
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Quadratic and cubic terms in the polynomial extrapolation represent small corrections of decreasing 
size. One can thus make some estimate of the gluon pole mass. A similar study on a lattice of double 
size, as was considered in Ref. fl^ j , would produce more than four times as many points in the same 
interval, and should allow for an accurate determination of the pole mass. 
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FIG. 7. Same as Fig.^ for the Laplacian gauge. The reduced vertical scatter of the data at a given momentum 
indicates a superior restoration of rotational symmetry. The reduced curvature as a function of improves the 
stability of the pole with respect to the type of extrapolation. 
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FIG. 8. Time-slice gluon correlator, in Laplacian gauge, at /3 = 5.8 and 6.0. The dashed lines show Corn- 
wall's model fitted to Dijf') after the cylindrical momentum cut; the solid lines are direct fits to the time-slice 
correlators, excluding the first few time-slices. 



We also measure the correlator of the gluon field averaged over a time-slice. Namely, we measure 



a— 1 /i 



(4.21) 



which is displayed in Fig.^. At large time separations this correlator should decay exponentially like 
exp{—mpoiet), giving us another approach to extracting the pole mass. We use this observable to perform 
a crosscheck on this mass, and as a further study of the systematic errors in its determination. This 
correlator is measured on the same configurations as D{q^), so it contains no additional information. But 
the same information is given a different weight, so that a fit to C{t) will give different results than a fit 
to D~^{q'^), especially after the cylindrical momentum cut. Therefore, we fit Cornwall's model directly 
to C{t) instead of D~^{q^). Remarkably, the difference is rather small, which attests again to the 
soundness of the model. The dashed lines in Fig.||show the original fit of Cornwall's ansatz to D^^{q^), 
which already provides a fair description of the data. The solid lines represent a direct fit of the same 
3-parameter ansatz to C{t), excluding the first few time-slices which otherwise completely dominate the 
fit. The fit started from t — A and i = 2 at /? = 6.0 and 5.8 respectively, which amounts to discarding 
similar intervals in physical units. Given the 3 fitted parameters, one can then solve D~^{q'^) = 



numerically, with D{q^) as per eg. ( 4. 17 ). The corresponding pole mass varies little from one fit to the 
other, and remains roughly constant in physical units aX (3 — 5.8 and 6.0. Also, a model-independent 
extraction of the pole mass, by measuring the effective mass TOe//(i) = —Ln{C{t + 1)/C{t)), gives a 
consistent value. Taking these results into account, together with the quadratic and cubic extrapolations 
displayed in Fig.^, we estimate the pole mass to lie in the interval [500, 785] MeV, where we used 
a-i(/3 = 6.0) = 1.943 GeV to convert to physical units ||6| with ^ = 440 MeV. The lower bound is 
given by the infrared mass, which corresponds to a linear extrapolation of D~^((f')] the upper bound is 
provided by the largest value obtained when fitting to our data Cornwall's model. A reasonable central 
value is 640 MeV, which corresponds to Cornwall's extrapolation in FigjTj 
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We have performed a similar exercise for the Landau gauge. The fit of Cornwall's model to D{q^) or 
C{t) is quite satisfactory, but the equation D~^{q^) — gives a complex pole, far from the real axis. 
Note that Ref. |3^] also finds oscillatory behaviour for the time-slice correlator in 3d SU{2) theory fixed 
to Landau gauge, reflecting a complex pole. This disagreement with the Laplacian gauge is puzzling, 
since one expects the pole to be gauge invariant. Possible causes include the inadequacy of the Landau 
gauge fixing procedure on the lattice, or finite-size effects. Larger volume studies, currently under way, 
should elucidate this issue. 

5. CONCLUSIONS 

We have evaluated the gluon propagator using the Laplacian gauge which avoids lattice Gribov copies. 
We extracted the transverse part of the gluon propagator and verified its scaling in this gauge. Examining 
the scaling and volume dependence of the zero- momentum propagator I'(O), we reached the conclusion 
that it is a constant beyond a lattice size of ~ 0.8 fm. This size is consistent with the characteristic 
length scale determined from I?(0) itself as the range beyond which the gluon field decorrelates in this 
gauge. 

Among the various proposals for the transverse propagator which arc physically founded, Cornwall's 
model [|o| provides a reasonable fit to the lattice results over the whole momentum range. We find it 
satisfying that the lattice data seem to favour a model with a dynamically generated mass. 

By looking at the inverse propagator D~^{q^) at small momenta, we see that the Laplacian gauge 
is superior to the Landau gauge in its restoration of Lorentz symmetry on the lattice. Furthermore, 
it turns out that the inverse propagator is almost linear in in the Laplacian gauge. This allows 
for a more reliable extrapolation to < 0, as compared to the Landau gauge. We test a variety of 
extrapolation ansatze. They consistently yield a pole mass at 640 ± 140 MeV. 
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